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SOLUTIONS 
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1. Introduction 

In this paper, we continue to consider the 2- dimensional (open) Toda 
system (Toda lattice) for SU(N + 1) 

N 

(1.1) -A Ui = J2 a v eU "> 

for i = 1,2, ■•• ,N, where K = (aij)N X N is the Cartan matrix for 
SU(N + 1) given by 
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As a very natural generalization of the Liouville equation 
(1.2) -Am = 2e u , 

system (1.1) is completely integrable, which is well-known in integrable 
systems theory. The Liouville equation and the Toda system arise in 
many physical models. In Chern-Simons theories, the Liouville equa- 
tion is closely related to Abelian models, while the Toda system is 
related to non- Abelian models. See for instance the books [10] and [24] 
and the references therein. 

Though the Liouville equation has been extensively studied after Li- 
ouville [16], the precise behavior of convergence of its solutions has been 
rather well understood only in the last two decades, see for example 
[3, 4, 5, 6, 7, 8, 14, 15, 17]. Such a delicate study leads to many applica- 
tions in the Abelian Chern-Simons theories and mean field equations. 

To well understand the non- Abelian Chern-Simons model, we have to 
study the analytic aspects of the Toda system. It is natural to ask if we 
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can generalize delicate analytic results for the Liouville equation to the 
Toda system. However, the analysis of the Toda system becomes more 
difficult, because the basic analytic tool, the maximum principle, does 
not work. In [11], we established a Moser-Trudinger type inequality 
for the Toda system, while a rough bubbling behavior of solutions to 
system (1.1) was considered. See also [13]. Its bubbles -entire solutions- 
were classified in [12]. In this paper, we will prove existence results 
for solutions to the Toda system in various cases by using methods 
developed in [8, 17, 6, 20, 19, 14]. And we will give a much more 
precise bubbling behavior of solutions, which will be very useful in our 
further study of the Toda system. 

Let £ be a Riemann surface with Gaussian curvature K. We consider 
the following system 

(1.3) -A Ui = !>;,<>,, [-7-^ 

for the coefficient matrix A = (cLij) nxn-i the Cartan matrix of SU(N + 
1) and p = (pi,p 2 ,--- , Pn) with p^ > (i = 1,2, ••• , N) given con- 
stants. Here hi : £ — > K. is a given positive C 1 function for % = 
1,2- •• , N. System (1.3) is the Euler-Lagrange system of the func- 
tional 

I N r N r N r 

(1.4) J p {u) = - / a ij VuiV Uj + J2 ~ J2 P^ ^ / 

i,j=l J j=l J j=l Jj] 

in H := (/J 1 (S)) iV where (a u ) is the inverse matrix of A. In [11], we 
proved that J p has a lower bound in H if and only if pi < An for any i, 
which is the Moser-Trudinger inequality for the Toda system. From this 
inequality, we know that if pi < 4n (Vj) then J p iscoercive condition 
and hence J p has a minimizer, which certainly satisfies system (1.3). 
When one of the p^s equals 4n, the existence problem becomes subtler. 
For simplicity of notation, we consider only the case N = 2. 
Our first result in this case is 

Theorem 1.1. Let E be a Riemann surface with Gaussian curvature 
K and N = 2. And let p\ = 4n and p 2 e (0, 4ir). Suppose that 



m 



E 1 < i < N, 



(1.5) A log tn(x) + (8tt - p 2 ) - 2K(x) > for iGE. 
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Then J p has a minimizer u = (u\,u 2 ) satisfying 
(1.6) 



. h ie Ul \ ( h 2 e U2 
A Ul = 2pi ^4 1 - P2 iH 1 



-Au 2 = 2p 2 . | - - 1 - Pl 1 - 1 



/or pi = An and p 2 £ (0, 47r). 

Theorem 1.1 is obtained by using a result in [6], which is a refined 
result of [8] and [17]. For p 1 = p 2 = An, we also have a sufficient 
condition under which (1.6) has a solution. See Theorem 5.2 below. 

For the general case, we have the following compactness of the solu- 
tion space. 

Theorem 1.2. For any compact set Ai x A 2 C R+ x R +7 if there are 
two positive integers mi and m 2 such that Aj C (Anrrii, An {mi + 1)) 
for i — 1,2, t/ien £/ie solution space of (1.6) for (pi, p 2 ) G Ai x A 2 is 
compact. 

The proof of Theorem 1.2 follows from the study of the convergence 
of the sequence of solutions, which was initiated in [11] for the Toda 
system. Let x be an blow-up point of the sequence u k , i.e., there exists 
a sequence — > x such that m&x{u k (xk),u k (xk)} — > oo as — > oo. 
Define 

<7j = lim lim / e Ui . 

We call (01,0-2) a blow-up value at x. First it is easy to check that 
°~i + °"2 > 0. A local Pohozaev argument gives us a relation between 
o~\ and a 2 . 

°\ + °2 ~ °"i°"2 = 47r(cri + a 2 ). 
See [13] and [11] for the proof. In Proposition 2.5 below, we show 
that (pi,p 2 ) can only be one of (Air, 0), (0, Air), (An, 8ir), (8n, An) and 
(8n, 8n). It is clear that Theorem 1.2 is a direct consequence of Proposi- 
tion 2.5. In the first two cases, one of the u\ does not blow-up, another 
bubbles more or less like solutions of the Liouville equation, which has 
been extensively studied in the last decade. The last case contains a 
typical blow-up phenomenon of solutions of the Toda system. It is one 
of our main results to obtain a precise description of the bubbling be- 
havior for this case. To describe it, we assume that there is a sequence 
of solutions u k = (u\,u k ) of 

1 -Au k 2 = 2h k e u * - h k e u * 
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Here B r denotes a disk of radius r and center and h k converges in 
C 1 to a positive C 1 function hi for % — 1,2. Assume without loss of 
generality that hi(0) = h 2 (0) = 1. Suppose that u k bubbles off, i.e., 
max xeB2 {u k , u k } — > oo as fc — > oo. More precisely we assume that 

(1) is the only blow-up point of u k . 

(2) max aB2 -u^ — min 9B2 -u t fc < c for i — 1, 2. 

(3) J" B e u * < c for i — 1, 2 and any k. 

Assume that \ k = \ k := max fl2 u\ > max B2 u\ =: X k . Let X& be the 
maximum point of u\. Set v k = (v k ,v k ) by 



v 4 (x 



' ' - Ui(e k x + x k ) - \ k for i = 1, 2, 



where e fc = e 2 A fc. Clearly t> k satisfies in {x e M 2 | e^x + x fe G -B2} 

-Auf = 2^(e fe x + x fc )e^ - h k (e k x + x k )e v " , 
-Auf = 2h\{e k x + x fc )e"2 - /if (e fc a: + x fe )e^, 

and i>f < 0, v\ (0) = 0. Since we only consider the case that the bubble 
is a solution of the Toda system, we may further assume that 

(4) v\ (0) is bounded from below. 

Then, there exists a solution v° = (^1,^2) °f the Toda system 

-Aw° = 2e"? - e v °, 
-Av° = 2e v ° - e v °, 



L(» 2 ) 



such that v k — v® converges to zero in C\ 

From the classification result for entire solutions of the Toda system 
[12], which is a generalization of the classification result for the entire 
solutions of the Liouville equation obtained by Chen and Li [5], v° is 
obtained from a rational curve from S 2 to CP 2 . In particular, we have 

e* = [ e «S = 8vr. 

Now we state our main theorem. 

Theorem 1.3. Let u k = (u k : u k ) be a sequence of solutions to (1.7). 
Suppose that h k converges in C l to a positive function hi with hi(0) = 1 
fori = 1, 2 and that (l)-(4) hold. Then there exist two constants r > 
and c > independent of k, such that 

(1.9) \u k (x) -X k - v?(e?(x - x k ))\ < c in B ro 

fori = 1,2. 
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Theorem 1.3 gives a precise asymptotic behavior of a blow-up se- 
quence of solutions. When N — 1, Theorem 1.3 was proved by Y.Y 
Li [14] by using the method of moving planes and a symmetry of the 
entire solutions of the Liouville equation. As mentioned above, the 
method of moving planes, which is based on the maximum principle, 
does not work for the Toda system. And also, the symmetry used in 
[14] is not available for the entire solutions of the Toda system in gen- 
eral. Thus the method in [14] could not provide a proof for the Toda 
system. Instead, here, together with the geometry related to the Toda 
system, we will employ a more delicate analysis. This is a combination 
of methods given in [2] and [12]. 

Our last result is an existence result for a supercritical case with 
respect to the Moser-Trudinger inequality established in [12]. 

Theorem 1.4. Let S be a compact Riemann surface of genus greater 
than 0. Then for any p = (pi, p 2 ) with Pi G (0,47r) U (4n,8n) system 
(1.6) has a solution. 

The method to prove Theorem 1.4 follows [20]. See also [21, 19, 8] 
and also [13]. One can also prove the Theorem by computing a topo- 
logical degree of system (1.6). In fact, the results -especially Theorem 
1.3- established in this paper will be used in computing the topological 
degree for the Toda system. We will pursue this subject in forthcoming 
papers. 

In Section 2, we first recall basic facts about the convergence of solu- 
tions to the Toda system and we then show that the possible blow-up 
values of the Toda system are isolated. We recall the geometric inter- 
pretations of solutions to the Toda system in Section 3. Such a solution 
corresponds to a flat connection and its singularity to the holonomy of 
its corresponding flat connection. We will give a more precise bubbling 
behavior, Theorem 1.3, in Section 4. The proof is a fine combination 
of arguments presented in [2] and [12]. In Section 5, we prove Theorem 
1.1 by applying a result in [6]. In Section 6, we give an existence result 
for (1.1) with the Dirichlet boundary for a supercritical case, based on a 
method of Struwe [20] and the Moser-Trudinger inequality established 
in [11]. We give an asymptotic behavior of singularities of solutions in 
Appendix. 

Acknowledgement. A part of the paper was carried out while the third 
author was visiting the Center of Theoretical Science in Taiwan (CTS). 
He would like to thank the CTS for warm hospitality. 
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2. Bubbling behaviors 

The convergence of solutions of Toda system was studied in [11]. 
Here for the convenience of the reader, we first recall some basic facts. 

Proposition 2.1. Let Q be a bounded smooth domain in 1R 2 and u k = 

{u k ,u 2 ) be a sequence of solutions of the following system 

, . ( -An* = 2h\e u * - h 2 e u %, on Q, 

\ -Au k 2 = 2h k .e u * -h k e u *, on tt, 

with 

(2.2) f e< < C, and [ e u " < C. 
Jn Jn 

Set 

(2.3) Sj = {x G E|3 a sequence y k — > x such that u k (y k ) — > +oo}. 

Then, one of the following possibilities happens: (after taking subse- 
quences) 

(1) u k is bounded in L% c (Vt) x L% c (fl). 

(2) For some j G {1, 2} 7 u k in L^ c (Q), but u k — > — oo uniformly on 
any compact subset of Vt for j 7^ i. 

(3) For some i G {1,2}, Si 7^ 0, but Sj = 0, for j ^ i. In this 
case, u k — > — oo on any compact subset ofQ\Si, and either, u k 
is bounded in L^ C (Q), or u k — > — oo on any compact subset of 

n. 

(4) Si 7^ and S 2 7^ 0. Moreover, u k is either bounded or — > — oo 
on any compact subset of VL\{S\ U E 2 ) /or j = 1,2. 

Proof. The proof is given in [11]. The idea follows closely [3]. □ 

Remark 2.2. One can prove that when S\ 7^ 0, S2 7^ and Si\S2 7^ 0, 
n^ — > —00 uniformly in any compact set of Q\{Si U S2). See [3, 11]. 
However, the proof of Proposition 2. 5 below implies that when Si 7^ 
and S 2 7^ 0, &o£/i u k — > —00 uniformly in any compact set of Q\{Si U 
S2). 77ms is an improvement of (4). 

In this paper, we do not distinguish convergence and subconvergence. 
We may assume that there exist two nonnegative bounded measures 
/ii and /i 2 such that 

e u ^ — > J ipd^i as k — > 00, 

for every smooth function ■?/> with support in f2 and i — 1,2. A point 
rr G is called a ^-regular point with respect to /ij if there is a function 
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ip G C C (Q) , < ijj < 1, with ■?/> = 1 in a neighborhood of x such that 

ipdfij < 7. 
We define 

^•(7) = {x efl\x is not a 7— regular point with respect to fij}. 

One can show ^1(7) and ^2(7) are finite. And ^(7) is independent of 
7 for small 7 < 2n, see [11]. Furthermore, we have 

(2.4) Si = ^(7), for 7 < 2n. 

Proposition 2.1 implies that the blow-up points are isolated. Let 
be an isolated blow-up point of the sequence u k , i.e., there exists 
a sequence Xk — > such that maxjMj^j^^Xfc)} — > 00 as — > 00. 
Define 

<7j = lim lim / e Ui . 

r^O k— >oo 

We call (<7i, <j 2 ) a blow-up value. If one of u\ and a 2 is zero, then another 
is An. This can also be obtained from the following local Pohozaev 
identity. 

Lemma 2.3. We have 

°"l + °2 _ a i a 2 = 47r((7i + cr 2 ) 

Proof. See the proof in [11] and [13]. □ 

Hence we assume that <7j > for i — 1, 2. From (2.4) we know that 
<7j > 27T for i = 1, 2. This can be improved in the following 

Lemma 2.4. We have o~i > An for i = 1,2. 

Proof. A proof of this Lemma is given in [13]. Here for convenience of 
the reader, we give a proof. Assume by contradiction that <j\ < An. 
From Lemma 2.3, we have <r 2 < 8n. Let Xk be the maximum point of 
u\ and \k = maxwf the maximum value of u\. Since o~\ > 0, we know 
that Afc^ooas/c-^oo. Consider 

u k (x) =u k l {e- Xk/2 x + x k ) - X k , 

in il k '■= {x E M 2 \ e~ Xk l 2 x + Xk G £v„} for a fixed r > 0. It is clear 
that (u\,u k ) satisfies (2.1) in Vtk and Vtk converges to M 2 as k — > 00. 
Now we have three possibilities: 

(i) M2 blows up. Namely there is a point y G R 2 and a sequence y& 
such that yk y and u^iyk) — > +00 as — > 00. 

(ii) M2 uniformly converges to —00 in any compact set of R 2 . 
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(iii) u\ is bounded from above and there is a point x such that 
u\{x) > —C for some constant C > independent of k. 

We first consider case (i). We define S 2 to be the blowup set of u\. 
It is clear that S 2 is finite. In fact, one can show that the number of 
S 2 is I. 

Since u\{x) < 0, by using the equation for vJ^ alone, we can apply a 
result in [13] to conclude that for each blow-up point y G S 2 , the local 
mass around y would converge to An. Thus 

An#S 2 < / e**dx <a 2 < 8n. 
Jn k 

Hence #S 2 — 1. 

Suppose S 2 = 0. Then both u\ and u 2 converge to v\ and v 2 , which 
are a solution of the Toda system (1.1) with N = 2. But the mass for- 
mulas (1.8) implies a 2 > 8n, which is a contradiction to the assumption 
a 2 < 8n. 

Now suppose S 2 = {p}. Then it is easy to see u\ converges to u in 
M 2 \{p}, and u satisfies 

-Au = 2e u - An5 p . 
Let v — u — 2 log \x — p\. It is clear that v satisfies 

-Av = 2\x-p\ 2 e v . 
A result given in [23] (see also [2]) gives that 

2 / e u = I \x\ 2 e v = 16tt, 



JR 2 JR 2 

which implies that o\ > 8tt. This contradicts o\ < An. 

In case (ii), one can show that u\ converges in Cf oc (IR 2 ) to a u, which 
is an entire solution of (1.2) with finite energy, the classification result 
of Chen-Li in [5] gives us 

lim / e"i = lim / e"* > Air, 



which contradicts o\ < An. In case (iii), we can show that (ui,u 2 ) 
converges to an entire solution of the Toda system (1.1) with finite 
energy. Now the classification result for the Toda system [11] gives us 
that 



lim / e"i = lim / e u i > 8tt, 
again, a contradiction. □ 
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From Lemma 2.3, the dimension of the set of possible blow-up values 
is less than or equal to one. We will show that the possible values of 
°2) are in fact isolated. It might be not difficult to see that (Air, 0) 
and (0, 47r) are possible values. And (8tt, 8n) is also a possible blow-up 
value. On M 2 , the solution space of the Toda system is noncompact 
and the blow-up value is (8ir, 8ir). Now we have 

Proposition 2.5. The blow-up value of the Toda system (2.1) can only 
be one of (An, 0), (0,47r) 7 (An, 8n), (8n,An) and (8n, 8tt). 

Proof. We only need to exclude the case that one of o~\ and o~ 2 is greater 
that 8n. Assume by contradiction that a 2 > 8n. In view of Lemma 
2.3, we have a 2 < V2n and An < U\ < 8n. 
Choose r k such that 

(2.5) / e< = 8n. 

Since a 2 > 8n, it is easy to check that r k — > as k — > oo. Consider 
u!?(x) =u k i {2r k x) -21og2r fe , 

in B r -i. Let Si be the blow-up set of u\. 

If S 2 — and Si ^ 0, then -u^ converges in Cf oc (lR 2 \ 1 Si) to a solution 

of 

(2.6) -At; = 2e V0 - An <J p . 

It is clear that the number of Si is 1. Otherwise, <7i > 87r, which is a 
contradiction to our assumption. Hence as above we have 

/ e V0 = 8tt. 
Jr 2 

However by (2.5), 8n = f R2 e v ° > f B e v ° > lim f B e u * = 8n, a contra- 
diction. Thus Si = (j). If S 2 = Si= 0, then as discussed in the proof of 
Lemma 2.4, either u\ converges to an entire solution v of the Liouville 
equation and u\ converges uniformly to — oo in any compact domain of 
1R 2 , or (u\, u 2 ) converges to an entire solution of the Toda system. The 
latter case is clearly a contradiction to o\ < 8ir and the former leads 
to 

An = [ e d0 > [ e s " > lim I & = 8n, 

JR2 J B2 k^oo J B _ ± 
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a contradiction again. Hence S 2 7^ 0- By rescaling a factor of close to 
1 if necessary, we may assume that {1 e I 2 \x\ = 1} fl (Si U S 2 ) = 0, 

(2.7) u k (x) -> -00, locally in R 2 \(Si U S 2 ) as k -> 00, 
for i — 1,2 and 

(2.8) / e n " _> 8tt as fc -> 00. 
•UW<i} 

Here (2.7) holds by Remark 2.2. Now we use the Kelvin transformation 
to define 

v k (x) = u k i^^j -41og|x|. 

Since the Toda system is conformally invariant, v k = (v k ,v k ) satisfies 
the Toda system in M. 2 \B rk . We consider v k on := Bi\B Tk . The 
boundary of Q k consists of two components 

diVl k = {\x\ = 1} and d 2 Q k = {\x\ = r k }. 

On diflk, v k — > —00 as k — > 00. Let fi k = max^ fe v k . Since a 2 — 8n > 0, 
it is clear that fi k —> 00 as k — > 00. Let y k G fife such that v 2 (y k ) = \i k . 
We claim that 

(2.9) ^00 as — > 00. 

Suppose that \yk\/r k < d uniformly for some constant d > 1. We 
consider -u fe in -B r -i. We first have by integrating by parts, for any 

% 

for some constant e > 0. Here 'uf(r) = ^ f\ x \ =r u k d<r is the average of 

£tf over {\x\ = r}. 

Let /(r) = u k (r) +41ogr. From above we have f'(r) < — e /r. 
Hence, we have for r > 1 

(2.11) fi*(r) + 41ogr < f(l) = =: c fc . 

From above, we know c k — > — 00 as /c — > 00. Now by the definition of 
t/fc, we have for any |x| > 1 



u k (x) + 4 log |x| = v k (-^r) < v k (y, 



(2.12) 



4(r%) ~ 41og|yfe|. 
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Recall that there is a constant c d > 0, independent of k, such that 

(2.13) \u k (x)-u k (y)\<c d , 
for \x\ > d~ l , \y\ > d' 1 . Thus 

(2.14) \u k {x)-u k {y)\ <c d 

for \x\ > r k x d~ x and \y\ > r^d^ 1 . (2.14), together with y^j > d -1 , 
implies that 

(2-15) \u k 2 (T%) ~ u k 2 (j%)\ < c d . 

mr \yk\ 

(2.15) , together with (2.11) and (2.12), implies that for \x\ > 1, 

u k (x) + 4\og\x\ < u k (-%) - 4\og\y k \ 

< ^(f^) -41og|y fc |+ Cci 

ml 

< Cfc + Cd — > — oo as k — > +oo. 

Therefore, we have 

as k — > oo, which implies that <r 2 = 87r, a contradiction. Thus, (2.9) is 
proved. 

Consider a new rescaled v k = (v k ,v 2 ) defined by 
v k (x) = v k (e-» k/2 x + y h )-Hk 

in 

tt k = { x g R 2 1 e -^ k/2 x + y k e Q k }. 
From (2.9), we have two possibilities: 

(i) \y k \e^ k ^ 2 — > oo as A; — > oo. 

(ii) 1 2/A; I e A * fc / 2 is uniformly bounded. 

For case (i), we know that £l k — > R 2 . Arguing as in the proof of Lemma 
2.4, we can show that converges in a suitable topology to an entire 
solution to an equation like (2.6) or (v k , v 2 ) to an entire solution of the 
Toda system. In both cases, we have 

lim / e"2 > 47r, 



c k 







which, it turn, implies that a 2 > 12n, a contradiction. For case (ii), 
we assume that ykC^I 2 — > p G R 2 . In view of (2.9), we have that 
f2 fc converges to R 2 \{p}. In this case, one can show as above that v 2 
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converges to a function v in Cf oc (M?\({p} U Si)), where Si is the set 
of blow-up points of v\. One can check that v satisfies 

(2.16) -Av = 2e s ° - aA7r5 p - An ^ 5 *- 

Since v\ < 0, we have a > 0. It is clear that the number of Si is less 
than 2 because U\ < 8n. If Si = 0, together with a > 0, then we 
have J R2 e v ° > An, which implies that a 2 > 12n, a contradiction. Now 

assume that Si = {q}. We will show in the appendix that J R2 e v ° > An, 
a contradiction again. This completes the proof of the Proposition. □ 



3. Geometric interpretations 

Before we give a proof of the precise asymptotic behavior of solutions 
of the Toda system, we would like to recall the geometric interpretation 
of solutions of the Toda system (1.1). 

Let Q be a simply connected domain and u = (u±,U2,-- - ,un) a 
solution of (1.1) on Q. Define wo, wi, W2, ■ ■ ■ , wn by 

N 

(3.1) Ui = 2wi — 2wi-i for i £ I and ^^Wi = 0. 

It is clear that 

N 




w = — y^(N - i + l)ui. 

2(N+1) ^ K J 



i=i 

Now set Wi = Wi — i log 4 for i — 1, 2, • • • , N. It is easy to check that 



e 2(w!-w ) 



e 2(w\— w ) _|_ e 2{w2-w{) 



2{w N ) z2 = - e 2 ^"^-i). 
(3.2) is equivalent to an integrability condition 

(3.3) U- Z -V Z = [U,V] 
of the following two equations 

(3.4) 4>- 1 -<t>z=U 
and 



(3.5) 0' 1 • 2 = V, 
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where 



U 



and 



V = 



Set 



/( 



Wo) 







{w N )zJ 



+ 



\ 



( o 



3 UIl--ii) 



W N -W N -1 



W N -W N -1 



J 

\ 



(3.6) a = -(Udz + Vdz). 

Then, a is a one-form valued in su(N + 1) 
Frobenius Theorem, we obtain a map <fi : Q — 



With the help of the 
SU(N + 1) such that 



a 







As a connection on the trivial bundle Q x C N+1 
flat, i.e., a satisfies the Maurer-Cartan equation 

1 



Q, a (or d + a) is 



da H — [a, a] 
2 



0. 



When Q is not simply connected, we cannot apply the Frobenius 
theorem directly to obtain <fi. Let Q = B* = B\{0} be the punctured 
disk. We introduce the holonomy (see [18]) of an SU(N + 1) connection 
a on the bundle Q x C N+1 — > Q along around 0. Let (r, 9) be the 
polar coordinates. Decompose a as a = a r dr + a^dd. a r and ag are 
su(N + l)-valued. For any given r G (0, 1), the following initial value 
problem, 

^ + a e <p r = 0, r (O) = Jrf, 

has a unique solution r (6 ) ) G SU(N + 1). Here /<i is the identity 
matrix. If the connection d + a has no singularity at 0, the conjugacy 
class of r (2-7r) is trivial. If a is a flat connection on D*, then the 
conjugacy class of r (27r) is independent of r. In this case, denote the 
conjugacy class of r (27r) by h a . h a is called the holonomy of a. 

Proposition 3.1. Let u = (ui,u 2 ,-- - ,un) be a solution of (1.1) in 
R n \{0} with 

Ui(x) — fa log \x\ + 0(1), near 0, 
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where Hi > —2 for i e / = {1, 2, • • • , N}. Then the flat connection a 
defined by (3.6) has holonomy 



(3.7) 



h, 



/ e 2nifo ••• \ 

e 2nilSl 



\ e 27rift 
where f3 , Pi, • ■ • , (3^ are determined by 



J 



N 



(3.8) 



Pi - = --Hi (i e J) and ^ ^ = 0. 



3=0 



For the proof, see [12]. 

Now by using Proposition 3.1, we can generalize a result in [12]. 
Given an iV-tuple (/ii,/i2, )A*jv) with //, > —2 for any i and let 
■u = (ui, u 2 , • • • , mat) be a C 2 solution of the following system 

(3.9) 



iV 



- AUi = ^ a ij 



m 



3=1 



with 
(3.10) 

Let 



/ 

Jm. 2 



|d Ml e Ul < oo. 



rrii = — f |x| Mi e Ui and 7« = 

2tt J R 2 ^ 

The potential analysis gives that 



U.; 



(3.11) 



= — 7, log \x\ + a,i + 0(\x\ x ), near oo, and 



x 



Vu t = - 7i — + 0(|x|- 1 ) 
\xr 



near oo, 



for some constants and 

(3.12) 7i -^i>2, 

for all %. For the convenience of the reader, we will give the proofs of 
(3.11 and (3.12) in the appendix of the paper. 

Proposition 3.2. If pLi <0 for any i, then we have 

7; = 2(2 + //;), for any i. 
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Proof. Let u = (ui, U2, • • • ,un) with in = Ui + Hi log \z\. It is clear that 
u satisfies (1.1) in M 2 \{0}. Let a be the flat connection (3.6) obtained 
from the solution u. We have 

Ui(x) — /Xjlog \x\ + 0(1), near 

iii(x) = — (ji — fXi) log |x| + 0(1), near oo. 

Let Pi be determined by G /} using (3.8). By Proposition 3.1, 

the holonomy of a at is h a (0) given by (3.7). Now we compute the 
the holonomy of a at oo, h a (oo). To compute it, we use the Kelvin 
transformation to consider 

x 

Vi(x) = Ui(— ) - 41og|x|, iel. 

Clearly, v = (v 1 ,v 2 , ■ ■ ■ ,v N ) satisfies (1.1) on IR 2 /{0} and 

Vi = (7i - A*i - 4) log |x| + 0(1), near 0. 

Let a' be its corresponding flat connection. It is obvious that h a (oo) = 
h a '{0), which can be obtained again by (3.7) by replacing //$ by {7^ — 
fii — 4}, i.e, h a f(0) is a matrix of form (3.7) with decided by 

ft ~ PLi = -\(H - ^ - 4) and J2 ft = 0- 

Now the key fact is h a (0) = h a (oo), which implies that 

pi -P'i = l modZ, for i = 0, 1, 2, • • • , N. 

Hence, 7, — 2/Xj = 2rii for some rii G Z. By (2.11), we have n« > 1 
for any i. We now need a global Pohozaev identity, which is given in 
Proposition 3.3 below. By (3.13), we have 

^V 7i (2n, - 4) = 0. 

Since a' 3 ' > for any i,j (see [12]) and rii > 1, the previous formula 
implies that rii = 2 for any %. This proves the Proposition. □ 

Remark 3.1. To remove the condition < /or any i ; one may aave 
to relate the solutions of the Toda system to other geometric objects. 
For instance, a generalization of projective connections considered in 
[22] might be a good candidate. 

Proposition 3.3. 7, satisfies 

(3.13) ^o ij 7 i(7i-2(2 + ^)) = 0, 

where (a*- 7 ) zs iae inverse matrix of the Cartan matrix (a^-). 
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Proof. First for a given R > 0, multiplying (1.1) by x • Vu k and inte- 
grating over B R , we have a Pohozaev identity 

(3 14 \ ^ "'Sfl 1 J b r 

From (3.14), we have 



(3.15) 



" 2^° .L/fri dn + ^JJ Xl 



xr 3 e 3 x ■ nuj. 



It is easy to compute that 
(3.16) 



\Y, aM I x-V(V jUk )V jUi = - ( x-V(V jUk V jUi ) 

J B R ^ </B fl 

1 Jb r 

+jVa* / RVu k V Ui . 



Inserting (3.16) into (3.15), we get 

[ (2 + fij)\x^e^-J2 I R\ 

JBf, JdBu 



<£ I W g^J 



Vu k Vui) 



Now using (3.11) in the above formula and taking the limit R — > oo, 
we have proved the proposition . □ 

A geometric proof of the Proposition can be obtained by using an- 
other geometric object, projective connections with singularities, which 
was used by Troyanov [22] to classify conical metrics of constant cur- 
vature. One can check that a quadratic differential r\ = fdz 2 , where 
/ : M 2 /{0} -> C is given by 

N ■ 1 

/ = 5^ a 3k {{uk)zz ~ ^(uj) z ■ {u k ) z } 

j,k=i 
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is a projective connection with regular singularities at and oo. By 
counting and comparing the weights of the singularities, one can get 
the global Pohozaev identity. 

4. A PRECISE BUBBLING BEHAVIOR 

In this section, we prove Theorem 1.3. Let u k = (u k , u 2 ) be a solution 

of 

' -Au k = 1h\e< - h k e u " 
-Au k 2 = 2h k e u " - h k e 
satisfying conditions given in Theorem 1.3. Set 



(4-1) 1 " T * T u ^ B 2 



lim / e 



/ 



From Proposition 2.5, we have o\ = o 2 = 8tt. 

In order to show the idea of our proof, we first prove Theorem 1.3 in 
the case that h k = h 2 = 1. Then we will point out which steps should 
be modified for the general case. 

Proof of Theorem 1.3 in the case that h\ — h 2 — 1. We divide the 
proof into several steps. 

Step 1. From the above discussion, (1.9) is valid in B €kR for any fixed 
large number R > 0. Hence, we only need to prove that (1.9) is valid 
in (efc.R, r ) for some r > 0. By assumption (4), we have 

\u k (x k ) - \ k \ < c, 

for some constant c > independent of k. In the sequel, c will denote 
a positive constant independent of k, which may vary from line to line. 



Step 2. By the Green representation, we have for \x\ < r , 
(4.2) 

u k {x) - m k = I G(x, y) (2e u " (y) - e u ^ dy + 0(l) 

= — [ \og—^—(2e u ^ y) -e uk ^Ady + 0(l), 

where rh k = infQ B2 u k (x). In this paper, 0(1) denotes a term bounded 
by a constant independent of k. Hence, at Xk we have 

(4.3) \ k -m k = — [ log , 1 , (2e"ik> - e u ^A dy + 0(1). 

27T J Bro (x e ) Fe ~y\ V ' 
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(4.2) and (4.3) imply 

{AA)v k {x) = u k (x k + e k x) - u k (x k ) 

= — I log J^L ( 2e v *M - e ti(v)) dy + 0(1). 
2tt J b ,-i \y-x\\ ) 



Similarly, we have 



k 2 (x k +e k x)-u k (x k ) = i- jf ^togJ^L^W-e^) Wo(l). 



u 

- - 2- , 

r k 

Since \u k (x k ) — \ k \ < c, we have 



(4.5) v k (x) = — [ log -^L ( 2e v ^ - e v " iy) ) dy + 0(1). 
Set w;* = |(2ui + uf) and w k = ^(v k + 2uf). It is clear that 



* = J- / to gl -!^e^dy + 0(l) 
/ log^L e ^*/ + 0(l). 



wS= w ... M 



A standard potential analysis (see for instance [2]) gives 



(4.6) 



w k + ^\og\x\ 



< 5\og\x\+0(l), 



for a small 5 > and |x| large enough, where erf is the local mass 
defined by 



Thus, we have 

(4.7) \v k -m k \og\x\\ < 3(5 log |x| +0(1), 

where 

ra? = — and rrC, = — -. 

1 2tt 2 2tt 

As in [2], we choose a small constant 53 > such that for log(^) < 
\x\ < 1 



(4.8) |m*(x)+m*| = 0(l)(loge^)-\ 

where 

mf(x) = / e v ^ y) dy. 

J\y\<S 3 \x\ 
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This can be done, since from (4.7) we have 
(4.9) 

I e v ^ y) dy<c! \y\- m > +3S dy = 0(1) (log e^ 1 )^ 2 , 

by noting that m\ — > 4 as k — > oo. It is also easy to check that 



/ 



1 logy^-e^dy < OMQoge?)- 1 , 

\ y \>8 3 \x\ F-y| 

for \x\ > loge^ 1 . Therefore, we have 

= — [ log-^-e v ^dy + 0(l) 

2n J\ y \<s 3 \x\ \y-x\ 



= ±- [ \ og \y\ e »l(y) d y + 0{l) 
2rr J\v\<s 3 \x\ \y-A 



] \y\<h\x\ 
ra^log \x\ + 0(1), 



for loge t < \x\ < — . Hence, we have 
(4.10) \v*(x) +m 4 fc log| x\ \ < c, 

for loge^T 1 < \x\ < Similarly, we can show that 



(4.11) 



x 



Vff(x) + mf 



\x\ 



< c 



holds for loge^, 1 < \x\ < j-. 



By (4.10)-(4.11), we have 
(4.12) 



u*(x) = m\ log^ + (2 -mf) log - + 0(1) 

\%\ 



V«?(x) = -m* * +0(1) 

in B rn \B . -i. 

Step 3. We want to show that 

m * =47r + 0(l)(log-)- 2 

for large k. Here we will use the geometric properties of the Toda 
system, described in Section 3. 
Let wo,wi and W2 defined by 

2 

(4.13) 2(wi - Wi„i) = (i = 1,2) and = 0. 
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Let wi — wi — log 4 and w 2 = w 2 — 2 log 4. Since ui, u 2 satisfy the Toda 
system, we already know that 

(4.14) da + ]-[a,a} = 0, 



where a is defined by 
Here 



2' 

a := Udz + Vdz. 



{w ) z e wi - w " 
U = ( (w 1 ) z e W2 ~ w 
(w 2 ) 2 



and 



N* 

V = - ( e'" 1 -'"" (im) 2 
e W2 " Wl (w 2 )i 

It is clear that a is a one-form valued in su(3). Equation (4.14) means 
that the connection d + a is flat. Now we calculate the holonomy of 
d + a along 

{|a;| = (login 

where s > will be fixed later. Decompose a by 

a = a r dr + agd9. 

By definition, the holonomy h a of the connection a is the conjugacy 
class of g r (2n), where g r {6) is the unique solution of 

dg r 

(4.15) -de +ae9r = °' 

0r(O) = 0. 

From (4.12), it is easy to check that 

/ 

a = j biie~ ie -i(3\ b 2 ie~ ie 

(4.16) \ Me* e -i/3 2 fe 

+ 0(l)e fc log-J 
at |x| = (log^) s , where 

1 rrfi 

Since m 8 > 2 for all i, we choose s > such that s(mhijmf — 2) > 4. 



21 

Now we can compute the holonomy h a of a along {\x\ = (log j^) s } 

(4.17) e" 2 ™^ +0(l)(log-)- 2 /. 
V e~ 2m ^ J €k 

Since the holonomy h a must be trivial, we have 

ft = n t + 0(l)(logV 
for some integer rii, which implies 

(4.18) m k l =2n' l + 0{l){\og-)- 2 

for some integer n'^ From Lemma 2.3 we know — > 4 as — > oo. 
Hence, we have 

m* =4 + 0(l)(logi)- 2 , 

which, together with (4.10), implies 

v\ = -m k i \og\x\+ 0(1) 

(4.19) = -41og|x| + 0(l)(log — )" 2 log +0(1) 
= —4 log +0(1), 

for loge fc 1 < \x\ < j-. 

Step 4. We have proved that (4.19) holds in loge^ 1 < \x\ < j-. By 
Step 1, we have, for any large R > 0, v\ converges to i>° uniformly 
for \x\ < R. By a classification result for entire solutions to the Toda 
system [12], we have 

\v° + 4 log \x\\ < c 
for \x\ > R, where c > is a constant independent of R. Therefore, 

(4.20) \vf + 41og|x|| < 2c 

for \x\ = R and large k. Choose R large such that 

e v > < \x\- 3 



for \x\ > R and i — 1,2, and define 

w±(x) = —4 log \x\ ± (ci — ci|x| _ 2 ). 

It is clear that 

1 , ,- 5 
Aw± = +-ci|x| 2 . 
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By choosing R and C\ large, we have, from the maximum principle, 
that 

w-(x) < v k (x) < w+(x) 

for R < \x\ < e^ 1 . Now we complete the proof of Theorem 1.3 in the 
case h\ = h\ = \. □ 

Proof of Theorem 1.3 in the general case. Assume without loss of 
generality that h k (x) = hi(x) and hi(0) — 1 for % — 1, 2. It is clear that 
we only need to prove Step 3, i.e., 

(4.21) m k = 47r + 0(l)(log V 2 , 



where 



m? = — and mZ = 



2tt z 2tt 

When hi is not constant, (4.1) is not an integrable system. However, 
the same idea as in Step 3 still works. 

Define Wi and Wi as above. We define a connection a = —{Udz+Vdz) 
on the trivial bundle, where 



U 



and 

/ (w )- z 
V = - tfe w i- w ° ( Wl )- Z 

\ fk^-W! 

Here ff = h l J 2 {e k x + x k ). Note that now v k = (fi,ff) satisfies in 
tt k := {i G K 2 e k x + x k e B 2 } 





(4.22) 



-Av k = 2h 1 (e k x + x k )e v i - h 2 (e k x + x k )e v % , 
-Av k = 2h 2 (e k x + x k )e v " - h x (t k x + x k )e v \ 



Now the connection d + a is not flat. But it satisfies 

I , (/{%-e-™ 

F = da+-[a,a]= -(f?) z e wi - w ° (f^ze™ 2 '™ 1 I dzAdz. 



2' 



o -(n)ze W2 ~ wl o 



Here F = F r gdr A d6 is the curvature of the connection d + a. 

Let gv(0) be the solution of (4.15) for this connection d + a. We 
denote the conjugacy class of g r {2ii) by h r a . Now h r a may depend on r. 
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Computing as in Step 3, we have (4.16) and 

/ e- 2 ^o \ 

(4.23) h r a =[ e~ 2m ^ +0(l)(log-)- 2 /, 
V e- 2 **® ) €k 

at r = (log^) s . 

Now we consider gauge transformations. Let : M, 2 — > SU(N + 1). 
From the connection a, we can get a new connection by 

a = (p^dcj) + _1 o;0. 

For the new connection, we compute h r & . 

Lemma 4.1 ([18]). We have h r & = h r a . 

Using this Lemma, we can compute h r a by using a suitable equivalent 
connection. This new connection is chosen as follows. By solving in 

n e \B R 

d<t> 

— + a r (p = 
dr 

0(0,0) = /, 

we obtain a new connection, a = (p~ l d(f) + _1 o;0. Decompose the new 
connection as above a = dedO + a r dr. Clearly, we have a r = 0. Hence, 
we have 

(4,4)^ ^ f = 

where F = F r gdr A <i# is the curvature of the new connection. It is 
well-known that F = <p~ l F(j). 

Now we estimate F r6 . We claim that 



(4.25) e v " <C\x\-' s , for |a;| > i2, 

for some fixed R > The claim will be proved at the end of the proof. 
From (4.25), it is clear that 

\(fi) z e Wl - W0 \ < Ce k e^ < Ce k \x\- 3/2 , 

for \x\ < R. We have similar estimates for the other entries. Therefore, 
we have 

l^rel = \Fre\ < Ce k \x\- 1/2 for \x\ > R. 
On Br we have \F r g\ = \F r g\ < Ce k . Therefore, from (4.24) we have 

Mr, 0)| < Ce(\og- k ) s I, for r = (\og\) s . 
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Now it is clear to see that 

hl = 0(l)e(\og^) s I : 
which, together with (4.23) and Lemma 4.1, implies 



/t = n l + 0(l)(logi)- 2 , 



and hence 



1 



ml = 4 + 0(l)(log fc , ^ 

e K 

This finishes the proof of Step 3 and the Theorem for the general case. 

Now it remains to check Claim (4.25). The idea is similar to Step 4. 
We define for Id > R a function w 



w = — 3 log \x\ + (ci — ci|x| 



We have Aw = — jCi\x\~% . Since m\ — >■ 4 as k — * oo, we may assume 
m\ > 3 for large k. Applying the maximum principle, we get the 
claim. □ 

5. Existence: A critical case 

Let E be a Riemann surface with Gaussian curvature K. In this 
section and the next section, we will study the existence of solutions of 

(5.1) —Aui = ^PjQij j—^. - lj , infi l<i<N, 

for the coefficient matrix A = (a i3 -) nxn , the Cartan matrix of SU(N + 1) 
and p = (pi, p2, ■ ■ • , Pn) with p^ > (i — 1, 2, • • • , N) given constants. 
System (5.1) is the Euler-Lagrange system of the functional 

I N r N f N f 

(5.2) J p (u) = / a ij VuiV Uj + Pi u i - ^Pi lo g / V^'. 

i,j=i J j=i J j=i ^ 

in H := (if 1 (E)) 7V . When pj < An, it is a subcritical case with respect 
to the Moser-Trudinger inequality established in [11]. By the Moser- 
Trudinger inequality, the functional J p is coercive and hence has a 
minimizer. When pi < Att and some of the p^'s are equal to An, the 
functional J p still has a lower bound, but it does not satisfy the coer- 
cive condition. It is a critical case and the existence problem becomes 
more subtle. Even in the case that N — 1, the existence problem is 
a difficult problem, see [8, 17, 6]. Here we will use results in [6] to 
obtain a sufficient condition under which the functional J p achieves its 
minimum. For simplicity we only consider the case N = 2. 
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Proposition 5.1. For a fixed p 2 G (0, An), define 

X = {Solutions of (1.6) for p 1 G (0, An] and p 2 }. 
If (1.5) holds, then X is a compact set. 

Proof. If X is not compact, then we may assume that there is a se- 
quence {pk} with p k < An and p k — > An as k — > oo, and a sequence of 
solutions u k = (u k ,u k ) satisfying (1.6) with p 1 = p k and p 2 such that 

maxjmax «J (i) , max (x) } — » oo 

as — > oo. Since p 2 < An, by Proposition 2.5, we know that max«{ — > 
oo as k — > oo and M2 is uniformly bounded from above. Let u> fc be the 
unique solution of 

(5.3) = » (jgf - l) == /* 

with J w k = 0. Since -u^ * s uniformly bounded from above, w k is 
uniformly bounded in H 2,q for any q < 00. Hence u> fc converges to u>o 
in C 1,Q for some a > 0. 

Let = + u> fc and v k = u\ — 2w k . It is clear that v\ and ^ satisfy 

(5-4) -At,* = 2p k ( - 1 

and 

with /i* = h x e w . It is clear that v\ + 2v\ = 0. Set v k — v k - 
log J* s . It is clear that v\ satisfies 

(5.5) -Av k = 2p k (h k e*t - l) . 

In order to apply Theorem 1.1 in [6] to our situation, we need to show 
that h\ = h\e~ w converges to h\e~ w in C 2 ' a . Let p be the unique 
blow-up point of u\. As in [11], one can show that v k — > 8nG(p,-) 
and v k = -v k /2 -> -AnG(p, ■) in i? 1,9 (E) (g < 2) and in C? oc (E\{p}), 
where G(x, •) is defined by 

-AG(x, 0=^-1 

and J G — 0. We apply a result of [14] (see also [2]) to (5.5) and obtain 
that 

(5.6) \v k (x) - U k (x)\ < c, for |z -p| < r 
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for some constant c > and a small constant tq > 0. Here 
U k (x) = log 



e Xk 



and Afc = maxtij. Furthermore, one can prove that 

(5.7) |A fc -log J h k e v *\<c, 
and 

(5.8) \Vv k (x)-VU k (x)\<c 

for some constant c > and |x — p| < r . See [4]. From (5.6), (5.7), 
(5.8) and u\ = —\v\ + 2w k , we have 

|Ve^| <ce-5 Afc |Ve _ 3 B i fc | < Cl (l + e-^|Ve-^ fc |) < c 2 . 

Therefore, the elliptic estimates implies that u> fc converges in C 2 ' a to 
w° satisfying 

- 1 

One can check by applying the Pohozaev identity that the blow-up 
point p satisfies 

(5.10) Vlog/n(p) + VG(p) - Vw°(p) = 0, 

where G is the regular part of the Green function G. Note that 
Aw°(p) = p 2 . Since w k converges to w° strongly in C 2 ' a , we can apply 
Theorem 1.1 in [7] to (5.5) and obtain that 
(5.11) 

2 Pk - 8tt = 2h^ (p)e w °^ (A log h (p) + 8vr - p 2 - 2X(p)) A^e"** + o(e Xk ) , 

where A^ = max {if, which tends to oo. Since 2p k < 8ir, (5.11) leads to 
a contradiction to (1.5). Hence X is compact. □ 

Proof of Theorem 1.1. Choose a sequence {p k } with p k < 4n and 
p k — > 4n. From above, we know that for p k = (p fc , p 2 ) with p 2 < 4-7T the 
functional has a minimizer -u fc . By Proposition 5.1, we have that u k 
converges to u°, which is clearly a minimizer of J p with p = (4n, p 2 ). D 

Theorem 5.2. Let T, be a Riemann surface with Gaussian curvature 
K. Suppose that 

(5.12) min{Alog/ii(:r), Alog/i 2 (x)} + Att-2K(x) > 0. 

Then J p has a minimizer u = (wi,w 2 ) satisfying (1.6) for p = (pi,p 2 ) = 
(47T, 47r). 
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Proof. Let p\ < An be a sequence tending to 4tt. Applying Theorem 
1.1, we know that J p has a minimizer u k = {u\,u k ) satisfying (1.6) for 
p = (47r, p\). To prove the Theorem, we only need to show that u k 
converges strongly to u°. Assume by contradiction that u k blows up, 
namely 



maxjmax-u^, max-u^} — > oo. 



Let 



Si = {x G £ | 3 a sequence {xk} with lim Xk = x & lim u { (xk) = oo}. 

fe^oo fc^oo 

By Proposition 2.5, we know that S± f] S 2 = 0. Hence we have two 
possibilities: 

(1) Si = { p } and S 2 = or S 2 = {p} and Si = 0. 

(2) Si = {p} and 5*2 = {q} with p ^ q. 

For case (1), using the argument given in the proof of Proposition 5.1, 
we can compute 87r — p\ in terms of the blow-up speed as in (5.11) to 
get a contradiction to condition (5.12). For case (2), we consider a disk 
DonE with p G D and q G" D. u k = (u\, u k ) satisfies on D 

— Aui = 8it ( -J^- l)—p 2 (-j!-T- — 

-Au 2 = 2p 2 „ , - - 1 - 4tt ' 



On the boundary dD, we know that there are two sequences {a k } and 
{a k } with — > — oo for i = 1,2 such that — is bounded. Define 



-Aw fc = P2 



Je ^ 



with the Dirichlet boundary condition w k (x) = for x G 9D. Now 
consider v k = u k + w k — a k which satisfies 



/ h k e< 



where 



h\ = e <h\e- wk . 
Note that v^ dD is bounded. Define Vq by 



-Avq = 8tt in D 



Vn = 



) k on dD 
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and set v k = v k + Vq. It is clear that v k satisfies 



-At;* = Stt-^^ inD 
| s fc fc e«f 

t) fc = on dD 
with h k = h k e~ v o . It is also clear that 



h h * e 



as k — > oo. As in the proof of Proposition 5.1, we can show that 
converges strongly in C 2 ' a . Hence, again we can use the result in [6] 
for the Dirichlet boundary problem to get that the divergence of v k 
implies that 

(5. 14) lim (A log h - Av% - Aw k - 2K) = 

at the blow-up point p. Since we have —Aw k — > —An and —Avq = 8n, 
from (5.14) we have 

A log h x + 4vr - IK = 0, 
which contradicts (5.12). □ 



A direct consequence of Theorem 1.1 is 

Corollary 5.1. Let S be a closed surface of genus greater than 0. 
Suppose that E has constant curvature and hi and h 2 are constants. 
Then equation (1.6) has a solution for any p = (pi, p 2 ) with Pi < An. 



6. Existence: A supercritical case 

In this section we consider the case that p« G (0, 8n), but pi ^ An for 
any i. This is a supercritical case, in the sense that the functional (1.4) 
has no lower bound. The result obtained in this section is not related 
to the curvature of the underlying surface. Therefore, in this section 
we consider the existence in a bounded domain Q C M 2 . A similar 
result as obtained in Theorem 6.1 holds for a closed surface with genus 
greater that 0, which is Theorem 1.4. 

Theorem 6.1. Let Q be a smooth bounded domain and pi G (0, An) U 

(47r, 8n) (i = 1,2) two constants. Suppose that the boundary dfl of Q 
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has at least two components. For any positive function hi (i = 1,2), 
there exists a solution u = (ui,u 2 ) satisfying 

M U1 h 2 e u > 
Awi = 2p 1 - r - — p 2 - 



y ' h 2 e U2 /lie" 1 

Am 2 = 2p 2 Pi r i M1 , 

and Ui = on dfl (i = 1, 2). 

When N = 1, the result was proven in [9], see also [19]. Recently, 
Lucia and Nolasco [13] obtained a non-trivial solution of (6.1) under 
some conditions when hi is a constant for % = 1,2. 

Our method of the proof closely follows [9], which, in turn, is moti- 
vated by [20] (see also [21]). Hence, here we only sketch the main ideas 
of the proof. When one of the pi is An, the existence becomes more 
subtle. We will consider this case elsewhere. 

Proof of Theorem 6.1. When pi < An (i — 1, 2), the Theorem is true for 
any bounded domain. We first consider the case pi G (An, 8n) (i = 1, 2). 
We divide the proof into several steps. 

Step 1. We first define the center of mass of a function v G Hq(Q) by 

Assume, for the simplicity of the notation, that dQ = d + Q U cLfi has 
only two disjoint components. Define a family of functions 

(6.2) T : (-oo, +oo) -> Hq(Q) x H^(Q) 
with T(t) = (71 (t),^ 2 (t)) satisfying 

(6.3) J P (T(t)) -> ±00 as t -> -00 
and 

(6.4) vrtJcii) — >■ 9 ± Q as t — > ±00. 

The existence of such a family is guaranteed by pi > An (i = 1,2). 
Denote the set of all such families by X and define a minimax value 

(6.5) a p := inf sup J p (T(t)). 

rex t 

Step 2. The minimax value a p > —00. 

The proof of the step follows from an improved Moser-Trudinger 
inequality under a condition introduced by Aubin. 
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Lemma 6.2. Let S± and S 2 be two subsets offl satisfying dist(S±, S 2 ) > 
5 > and 70 G (0,1/2). For any e > 0, there exists a constant 
c = c(e, 5q, 70) > such that 

holds for all u = (ui,u 2 ) G -f^o(^) x Hq(Q) satisfying 

Is eUl Is eUt 

(6.6) 1 > 70 and r 2 - > 7o, 

fori = 1,2. 

Let 70 C Q be a closed curve in f2 enclosing the inner boundary of 
VL. Each curve starting from <9_f2 and ending at d + Vl intersects with 
7 . By Lemma 6.2, we can show that 

J p (u) > -c, 

for any u E Hq(Q) x Hq(Q) with center of mass m c (u) G 70 . See the 
argument in [9]. Hence we prove the step. 

Step 3. a p : (4tt, 8n) x (47T, 871") is non-increasing in the following sense: 
if p = (pi,p2) and p' = (p' l5 p' 2 ) with pi < p' x and p 2 = p' 2 , then 

«p > ay 

Pi ~ Pi ' 

This is easy to check. 
Step 4. Now fix P2 and define 

Ai = Ai(p 2 ) = {pi G (4tt, 8tt) I is differentiable at pA. 

Pi 

It is clear that Ai is dense, i.e., A x = [47r, 87r]. Following a method 
given by Struwe [20] and [21], see also [8], we can prove that for fixed 
p 2 and pi G Ai, a pltP2 is achieved by u = (ui,u 2 ), which is a solution 
of (6.1). 

Step 5. Now for any p = (pi,p 2 ) G (Air, 8tt) x (47r, 8tt), we have a 
sequence p k = (p^,p 2 ) with p^ G A x and -u fe = (u k , ii 2 ) satisfying (6.1) 
with p fc . 

Since p 2 7^ 47rm, by Proposition 2.5 we see that u k converges to 
u° = («5,m 2 ), which is a solution of (6.1). This finishes the proof for 
the case p« G (Air, 8n) (i = 1,2). 

When pi < An and p 2 G (47r, 8tt), F defined by (6.2)-(6.4) does not 
exist. We need to modify the definition. Consider T = (71,72) with 
(6.2) and 

(6.7) ra c (^ 2 ) — * d±VL as t — > ±00. 
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Then the same argument finishes the proof of the Theorem. □ 
Theorem 1.4 can be proven in a similar way. 

Furthermore, in this case we can compute the topological degree of 
the solution space as in [7]. First we define the Leray-Schauder degree 
for system (1.6). For the Leray-Schauder degree for the Liouville type 
equation, see [14] and [7]. Consider p = (p±, p 2 ) with ^ 0mod47r. 
Define an operator 



T p := A" 



2 

-1 




T p acts on H% x H^, where H% = {u G H 1 \ j u = 0}. Set X p the 
solution space of (1.6). From above, we know that X p is compact. 
Hence we can define the Leray-Schauder degree for (1.6) 

d p :=deg(I + T p ,B R ,0), 

where I is the identity and Br := {u = (ui,u 2 ) G Hq x Hq \ \\ui\\ + 
1 1^2 1 1 < R} for a large R > 0. It is clear that d p is well-defined. The 
homotopy invariance of the Leray-Schauder degree implies that d p is 
independent of h\ and h 2 - Furthermore Theorem 1.2 implies that d p = 
d p if there are two integers mi and m 2 such that pi G (47r(mj — 1), Aurrii) 
and pi G (47r(mj — l),47rmj). 

Theorem 6.3. Let S be a closed surface of genus g. Then we have for 
p = (pi,p 2 ) with pi G (0, Att) and p 2 G (47r, 87r) 

d p = 2g- 1. 

The proof follows from the argument given in the proof of Theorem 
1.1 and a result given in [7]. In a forthcoming paper, we will compute 
the topological degree for the Toda system for the general case. 

7. Appendix. Behavior of singularities of solutions with 

finite energy 

In this Appendix, we consider the asymptotic behavior of singulari- 
ties of solutions to the Toda system. 

Let u — (ui, u 2 , ■ ■ ■ ,%) be a solution of 

N 

(7.1) -Au i = ^2a ij e u i, for % = 1, 2, • • • , iV 

i=i 

in a punctured disk D* = D\{0} with 



(7-2) / e« 

Jd 



< oo. 
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By the potential analysis, one can check that there is a constant 7, 

(7.3) lim — — —. — - = 7j, 

M-o-log|a;| 

for each i. (7.2) and (7.3) imply that 7i < 2 for any i. Furthermore, 
one can check that 



(7.4) lim ' ' , = 0. 

\x\->0 log |X| 

Proposition 7.1. For any i, we have 

7* <2. 

Proof. Suppose that there is j £ I such that 7^ = 2. We claim that 
7j_i = 2 or 7j +1 = 2. Assume by contradiction that 7j_i < 2 and 
7j + i < 2. We assume that neither j = 1 nor j = A/". (For the cases 
j = 1 and j = A/", the argument is the same.) Consider in D* 

-Auj = -e**- 1 + 2e u J - e u '' +1 =: /. 

Since 7^ = 2, 7j_! < 2 and 7^+1 < 2, we know that f(x) > in a small 
punctured disk D* 6 = D s \{0}. Set 



v(x) = / log|x - 
27r Jd s 



y\f(y) 



and w = uj — v. It is clear that Av = — f and Aw = 0. One can check 
that 



which implies 



lim — — = 0, 

|x|— »o — log \x\ 



W(X) ,. Uj —v 

lim i-f— = lim —2 — — = 2. 

|x|-*0 — log |x| |x|-»0 — log |x| 



1 



Since w is harmonic in Dg, we have w = —2 log |x| + wq with a smooth 
harmonic function wq in Z)^. By definition, we know v > 0. Thus, we 
have 

e u * = / e^" > f j-^e wo = 00, 

'D 4 JDg Jd s 

a contradiction. A similar argument implies that if 7, = 7^+1 = • • • = 
jj + k = 2, then either 7j_i = 2 or jj+k+i = 2, by using the equation for 
\{uj + Uj + i + • — h Uj + k). Hence we can show 7« = 2 for all i. Now we 
consider u = j^(^2f =1 Ui)- It is clear that u satisfies 

-Aw = ^(e ul +e U2 ) > 
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with f D e w < oo and lim^^o _ lc ^ ^ = 2. The same argument as above 
gets a contradiction. Hence 7$ < 2 for any i. □ 

Proof of (3.11) and (3.12). Now we apply the Proposition to solutions 
of (3.9) with condition (3.10). Let u = (u\,U2,--- ,un) be such a 
solution. The potential analysis gives us 

1 N f 

(7.5) lim = - 7i = ~^J2 a v / \ x \ HeUj - 

We consider u = (Hi, u 2 , ■ ■ ■ , Un) with f D e Ui < 00 for any i and 

Ui = Ui(-^-r) - (4 + /ij) log \x\. 

\x\ 

It is easy to check that u satisfies (7.1) in D* with 

,. uAx) 

hm — - — — - = 4 + ^ - 7j. 

|x|->0 — log |x| 

Applying the Proposition, we have 

4 + fa - 7i < 2, 

for any i. This is (3.12). From the above inequalities and the potential 
analysis, we have 

Ui = -(4 + - 7i) log |x| + 0(1), 

from which we have the first formula in (3.11). The second formula 
follows from the potential analysis. □ 

Lemma 7.2. Let u e C 2 be a solution of 

(7.6) -Am = 2\x\ 2a \x - q\ 2 e u 
with a > and 

(7.7) !■■=—[ \x\ 2a \x - q\ 2 e u < 00. 

2vr Jr2 

t/ien 7 > 2 + a. 

Proof. We show this lemma following closely [23] . Set 

w = — [ (log \x-y\- log(|y| + l))2|y| 2c % - g| 2 rfy. 

It is easy to check by the potential analysis that u = w + c for some 
constant c and 

w > —27 log \x\ — Ci, 
for some constant C\ > 0. This, together with the finiteness condition 
(7.7), implies that 7 > 2 + a. □ 
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Remark 7.3. In fact, all solutions of (7.6) and (7.7) can be classified. 
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